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Abstract 

By methods of quasiclassical asymptotics the behaviour of the integrated density 
of states for ID periodic nanostructures at the zero bias limit is studied. It is 
shown that the density of states at the zero bias limit has no regular limit while 
the integrated density of states has. The rigorous proof of this phenomenon given 
in the paper is based on a novel approach for the quasiclassical asymptotics on the 
spectrum of the Stark- Wannier operators. A connection of this phenomenon with 
the zero bias limits of the current through the nanostructures and their conductivity 
is briefly discussed. 

1 Introduction 

Some experimental works with tunneling devices have observed an increasing con- 
ductivity near zero bias [1]. This "zero-bias anomaly" has been explained in a 
number of papers by different physical phenomena [2, 3, 4]. Some of these explana- 
tions are based on the electron density of states and its behaviour near zero bias. 
In the present paper we study the behaviour of the density of states and integrated 
density of states at zero bias and show that the first one has no regular limit while 
the latter has. 

Because the current through the device and, in turn, the conductivity are defined 
through the integral from the density of states and some smooth function, it means 
that we should look for a "weak limit" of the density of states. We proof rigorously 
that the integrated density of states tends to an appropriate limit as the strength of 
the applied electric field goes to zero. This prove is the main novelty of the paper. 
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The main mathematical objects in our constructions are the spectral projections 
kernels of the operator 

d 2 

which describes ID periodic structure embedded in an external electric field with 
the strength e > 1 . The function p is assumed to be a periodic one. The operator 
H has been an object of the intensive investigations (see [14, 15, 16, 17, 18, 19] 
and references therein). However to our best knowledge the behaviour of spectral 
projections and of the integrated density of states which is expressed in their terms 
have not been considered as e — > 0. 

In papers [16, 17, 18] there has been performed the quasiclassical study of the 
equation Hg = Eg. These results have been used in [14, 15] for the solution of a 
number of spectral problems related to the operator H. For instance, there has been 
studied the behaviour of the residues of the analytical continuation of the resolvent 
kernel of H and proved that in an appropriate sense the ladders of resonances are 
concentrated e^Oon the spectral band of the operator H$. However the methods 
used in the above papers do not allow to study the behaviour of the operator H 
characteristics on the spectrum which fills the whole axis. In the present paper we 
overcome the difficulties related to the asymptotics constructions at real energies 
and hence provide the mathematical basis for studying the integrated density of 
states and the conductivity near the zero bias limit. 

One of the advantages of the approach proposed in the present paper is that 
we explicitly take into account in the leading order of e the tunneling probabilities 
through spectral gaps for the device embedded into the electric field. These prob- 
abilities are known [17, 18, 19] to be exponentially small w.r.t. to e however they 
strongly influence the behaviour of the integrated density of states. 



2 Preliminaries 

In this section we recall some facts concerning the one-dimension Schrodinger oper- 
ator 

with a periodic potential p(x + a) = p(x) and the Hamiltonian H given by (1). 

The spectrum of H in L 2 (k) consists of bands [X21-2, ^21-1] separated by gaps 
(A2j, A2i-i) [11]. In the present paper we assume that the number of gaps in the 
spectrum is finite, i. e. I — 1, 2, . . . , TV. Let us note that the class of the TV- 
gap potentials has been described completely [12] in connection with the algebraic- 

1 It might be the ID periodic nanostructure generated by the stimulated Mott-Peierls transition 
on quantum wire (see [6-10] ) with the bias applied 
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geometrical integration of the Korteveg-de Vries equation. However in the present 
paper we do not need the knowledge of the explicit form of these potentials but only 
make use of the fact that number of gaps is finite. 

The integrated density of states for the operator H is given by the formula [5]: 

1 ° 

N(X) = - ( e(x,x,X)dx, (2) 
a J 
o 

where e(x,y, A) is the spectral function of H : 

m 

e(x,y,X) = Y,ej{x,y) 

if A belongs to the m + 1-gap (A2 m , A2 m +i) and 

1 ^(x,X)i^y~X)] 



e(x,y, A) = ^2ej(x,y) + - j 

J_ 1 Mm 



X{k) 



dX 



if A belongs to the (m + l)-th band [A 2m , A 2m +i]- Here ej(x,y) is the kernel of the 
spectral projection of H on the j-th band Aj = [A2J-2, A2J-1]: 



1 r mx ,XMy,X)] dx 

TV J X(k) 



The Bloch functions ip(x, A) here are assumed to be normalized by the condition 

1 a 

- U(x,X)i^X)dx = l. (4) 
a J 


By X(k) we denote ^A(fc) where X(k) is the dispersion function which give the 
dependence of the spectral parameter A on quasimomentum k [13]. Recall that the 
Bloch function ip(x, A) is a quasiperiodic one and ip(x + a, A) = e %kx i\){x, A). 

Thus the integrated density of states for the operator H is nondecreasing func- 
tion which equals at A < inf a(H ) = Ao and is constant on each gap Aj = 
(A2J-1, X2j)- It is known that 

N(X) = M 

7T 

where k(X) = ^tk(X) and k(X) is the global quasimomentum introduced in [13]. The 
quasimomentum k(X) is real in bands and has nontrivial imaginary part inside gaps. 
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We introduce the integrated density of states for the operator H given by (1) as 
follows 

1 a 

N e (E) = - [e £ (x,x,E)dx (5) 



where e e (x, y, E) is the corresponding spectral function. In [14, 15] it has been shown 
that e e (x,y,E) has the form 

1 J g(x,E')g(y,E') ^, 
eJx.y.E) — — / . . . rz dE . 6 

— oo 

Here g(x, E) is a solution of equation Hg = Eg which satisfies the condition g — > 
as x — > — oo. Let us normalize it by the asymptotic relation 

g {x, E) ~ c e -Ml^l 3/2 ( 7 ) 



where the real constant c will be fixed later. The function M(E) is the so-called 
Jost function and is defined by the equation 

MiE) = ,8) 

Here W[g,f] = (?/' — </'/ and f(x,E) is the solution of the same equation having 
the oscillating behaviour as x — > +oo and normalized by the asymptotic relation 

v ; v/2 v 7 ^ + ex 

From equations (6) and (8) it is clear that e £ (x,y,X) does not depend on the 
choice of the constant c in (7). However it is convenient to fix the normalization 
of the Jost function M(E) unambiguously. Since we shall study the behaviour of 
N e (E) end hence of e £ (x,x,E) as e — > it is enough to fix the constant c in the 
leading order of e. Let us set 



c = exp 

where 



/ A \ A 

- J \k{\)\d\- J V^\d\ \ - J $(A)dA 

Voo — oo / — oo 



(I + 0(e)) (9) 



$(A) = i(«9 A 0(x,A),0(x,A)), (10) 
a 

4>(x, A) is the periodic part of the Bloch function if) normalized by the condition (4) 
and (•, •) is the standard inner product in L 2 (0, a). Since \k(X)\ ~ \J\X\ as |A| — > oo, 
the difference of integrals in (9) converges. 

To study the behaviour of the integrated density of states for Bloch electrons in 
electric field (5) at the limit of vanishing strength e we have to study the asymptotics 
of the spectral density \M(E)\~ 2 and of the solution g(x, E) at real values of E. 
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3 An asymptotic behaviour of the spectral den- 
sity. 



In [14, 15] there has been obtained the asymptotic formula for M(E) which allowed 
to study the structure of roots of M(E) in the lower half-plane (resonances) and 
to estimate their imaginary parts. However the formulae obtained there can not be 
effectively used to study the behaviour of the spectral density \M(E)\~ 2 as e — > 
on the spectrum and hence the behaviour of the integrated density of states (5). 
As it has been mentioned in Introduction in the present paper we propose a novel 
method of analysis of the Jost function. The advantage of the asymptotic formulae 
obtained below is that although we study the spectral density only in the leading 
order of e we artificially preserve a number of terms which are exponentially small 
with respect to e. However it is these terms which allow us to effectively study the 
behaviour of the integrated density of states as e — > 0. 

In the case of iV-gap potential the Jost function M(E) is described in the leading 
order of e as follows 



M(E) = e*F N (E)[l + 
The function F N (E) is given by the matrix product 




(11) 




(12) 



where F = —i and £ is some phase factor. 
The matrices 




(13) 



I — 1, . . . , N, have the entries 




(14) 




(15) 



Here 



En 



(16) 



oo = 

ea 



is the rescaled energy and ui is the constant: 




■nl/a 



(17) 



n(l-l)/a 
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where &(k)dk = $(A)dA with $(A) 
formed over the /-th Brillouin zone 
band A; under the mapping A;(A). 
The constant Pi is given as follows 



given by (10). The integration in (17) is per- 

tt(Z-I) ttZ] 
a ' a 



which corresponds to the /-th spectral 



P l = e" 2 ^ (18) 

with 

IH = -J |9*(A)|dA. (19) 

A; 

Here A/ = (X21-1, X21) is the /-th gap in the spectrum of Hq. It should be noted that 
Pi is the probability of transmission of an electron under the influence of applied 
electric field from the conducting band A^ to the band A i+ i ( see [16, 17, 19] and 
references therein). 

The phase factor £ in (11) also can be described explicitly. However we do not 
need its knowledge since to study the integrated density of states N e (E) (5) we need 
\M{E)\ only. 

In what follows we show how to simplify the matrix product (12). Let us intro- 
duce the parameter 

oi = iln(l-fl)- 1 (20) 
On the basis of (18) and (19) one sees that 

<n ~ i e - 27rw ^ as e -> (21) 

Thus the quantity 07 can be considered as an additional parameter of the problem 
which is exponentially small in e. Let 

a = max{o"/}. 

Then the following statement is valid. 

Proposition 1. Let the function Fi(E) be defined by the matrix product (12). 
Then in the leading order of a 

This estimate is uniform w. r. t. E el. 

Proof. This statement can be proved by the mathematical induction method. 
For I = 1 one has 

Fi(£) _ e -^-i(^i) _ e i(«-wi) ^ 1 _ e 2i(«-«i) + Q( a j 

~F[{E) ~ e i(w ~ wi )~ 2cr — e^^^^i) ~~ e 2i ^~^ + 0(a) ~ 
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Here we used the fact that F = —i and in the second equality we expanded e~ 2cr 
into the Taylor series as a — > 0. 

Now suppose that the statement is valid for I — 1. Since (F h F t ) T = X/(F;_i, F/) T 
one has that 



Fi — i 



-Act, 



Hence 



Fi 



-1/2 



e i(w-wi)-2cr i jp _|_ e -i(a;-a)i) ^ 



( c~ 2a ' Fl ~ 1 + e 2i ( a, ~ an ) 

c 2i(u-un) F l-i _|_ J 



(23) 



(24) 



Now expanding e~ 2cri as 07 — > and using the induction supposition we obtain (22). □ 
On use of this proposition one can obtain the asymptotics of F\ as a — > 0. 
Lemma 1. On the real axis E the functions F t ,l = 1,...,N, in the leading 

order of a are given as follows 



where 



Fi(E) = e l ^ E) ft (l - 



i(E) 



-4<Ti 



-1/2 



e 2 *'^'-^ - l) [1 + 0(a)], (25) 



z i=i 



(26) 



Proof. Take out from square brackets in (23) the multiplier Fj_i, and then use 
the recursion procedure. This yields: 



1 -1/2 1 



3=1 3=1 



,J(uj—u)j)—2<jj I ^_J_ 



+ e 



-i(u>-u>j) 



Now the statement of the lemma is obtained easily on use of the asymptotic relation 
(22) for the ratio (Fj^/FjTi). □ 

Taking now the absolute value of the r.h.s. of (25) one easily obtains \F N \. To 
write down it in a compact form let us introduce the function 



, x tanhcr 

o (to, uo , a) = — s- 2 v, 7 

sin (uj — uo') + tanh a cos 2 (u — uo') 



and the functions 



En 

5j(E) = 8(—,u>j,<Tj), j = l,...,N 



(recall that ujj and Oj are given by (17) and (20) respectively). Then 

1 N 



\F N (E)\ 



(27) 
(28) 

(29) 
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Now on use of the asymptotics (11) one obtains the main statement of this section. 

Theorem 1. In the leading order of e the spectral density \M(E)\~ 2 is given 
as follows 



N 



J'=l 

where the functions 8j are defined by (28). 



\MlEW-U^( E )^ + 0(e)), (30) 



4 Weak limit of the spectral density. 

The asymptotic formula (30) obtained in the previous section allows to prove that 
the spectral density \M(E)\~ 2 in a weak sense as e — > tends to unit. This important 
fact will be essentially used in the study of the integrated density of states (5). 

Theorem 2. Let AE be an interval of the real axis E and f(E) be an arbitrary 
smooth function defined on it. Then 

H [fgg* = [new- (3D 



£^0 



\M(E)\- 



AE 1 x ' 1 AE 

On use of the asymptotic relation (30) it is enough to show that 

/ ( n f( E ) = J f( E ) dE - (32) 

AE V'- 1 / AE 

Will shall prove the latter fact by several steps. 

Proposition 2. Let Auj n = [uo 1 + Tin — ir/2, iv' + nn — tt/2], uj' be the fixed point, 
n be an arbitrary integer and f{uj) be a smooth function on Au n . Let 5(u,iv',cr) be 
defined by (27). Then 

lim — f 8(ou, u', a)f{uj)duj — f(ou' + im) (33) 

Au n 

Proof. Make the change of variables 

z = tg {uj — uj' — nn) 

and denote 

a = tanhcr, f(z) = f(z + uj' + nn) 

Then 

lim— f 8(u, u', a) f (u)dw — lim — f — dz. 

Aqj„ — oo 
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Due to the fact that cr/(cr 2 + z 2 ) as a — > is the 5-type sequence the r.h.s. of the 
latter equals /(0) = /(a;' + 7m). □ 

Lemma 2. Let AE be an real interval, f(E) be a smooth function and 5j(E) is 
given by (28). Then 

lim J 5 j (E)f(E)dE = J f(E)dE. (34) 

AE AE 

Proof. On use of the definition of Sj one has 

J 5 1 (E)f(E)dE= £ -^ J 8{u,Uj,aj)f{u)dw, (35) 

AE Au) 

where 

Au = and f(u>) = /( ). 

ea 7r 

Let us split the interval AE on subintervals which have the length ea and the 

middle points 

( n) = ea +£an =a f + £ $(^.)]^ + £an 

IX ix J 

A, 

(here n e Z is such that i?] n ^ G AE 1 ). Under this procedure the interval Acj = 
is splitted into subintervals Auj^ = [uoj + nn — ir/2, Uj + nn — tt/2] where Uj is 
given by (17). Thus the integral in the r. h. s. of (35) can be presented as the 
sum of integrals over Ato^p. According to proposition 2 each of them as aj — > 
tends to 7xf(ujj + nn) = nf(Ej + ean). Thus the r. h. s. of (35) takes the form 
ea J2' n f(Ej+ean) where summation is performed over such n that Ef ] = Ej + ean e 
AE. Obviously as e — > this sum gives the integral in the r. h. s. of (34). □ 

The statement of lemma 2 differs from that we need to prove, namely equation 
(32), by the change of the single multiplier Sj(E) by their product. In this connection 
we have to prove that in appropriate sense for some j 

N 

llS^E) ~ 6j(E) as e^O (36) 

1=0 

Due to the singular behaviour of 8j(E) as e — > the direct analytic proof of (36) 
meets serious difficulties. Therefore we shall use the relation (29), take into account 
that the function Fn(E) is defined by the matrix product (12) and perform more 
detailed analysis of the latter. 

Together with functions F N (E) let consider the functions F$(E) which we define 
as follows 

Fn Fff \ I F,_! if,., \ 

TV iff " X '\TU -iFu) 
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On use of (14), (15) one has that 



detX t 



Hence 

and therefore 



TP Z?(0 TP pW 



Fi-i 



-i\Fi-A 



-20 i ^ 



(37) 



Now one can easily obtain the asymptotic formula for F$ (E). To this end repeat 
the proof of lemma 1 taking into account the additional multiplier — | in definition 

of F$(E). The result is the following statement. 



Proposition 3. In the leading order of a the functions F$(E) have the form 



N I/O 

F^{E)=ie^^X\(l-e-^) ' (. 



,2i(u>-u)i)-2ai 



+ 



N 

J] (e*("-"i)-*°i - l) [l + 0(a)]. 



Compare now the obtained formula for F$ with the equation (25) for F N . Ob- 
viously 

Inserting this representation into (37) and performing simple calculations one ob- 
tains the following statement. 

Lemma 3. On the real axis E in the leading order of a 

2 

= S l (E)(l + 0(a)) (38) 



F N (E) 

when Si(E) is given by (28). 

The estimate (36) is just the simple consequence of (38). Indeed, set / = 1 and 
note that F = — i. Then 



JV 



n^(£)[i+o(<7)] 



F N (E) 



5 1 (E)[l + 0(a) 



(39) 



Here the first equality is due to (29) and the second follows from (38). 

Now the statement of the theorem formulated in the beginning of the present 
section is the direct consequence of lemma 2 and the relation (39). 
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5 Nonuniform asymptotics of the solution g(x, E) 



In this short section we essentially use the results of the paper [14]. However we 
need this section since the nonuniform w. r. t. E asymptotics of solution g(x,E) 
combined with the behaviour of M(E) studied in the previous section define the 
asymptotic behaviour of the integrated density of states (5). 

Quasiclassical analysis of equation Hf = Ef performed in [14, 15, 16, 17, 18] 
leads to an important link between the dispersion function \(k) of the operator H 
and the variable E + ex which naturally arises in the asymptotics of g(x, E) given 
by (7) as x — > — oo. This the link has form 

X(k) = E + ex (40) 

and is valid for all x. At any fixed x this link can be considered as identification 
between the real axis E and the spectral axis A of the operator H . Let us denote 
by 

Aj(e) = [\21-2 - ex, A 2 «_i - ex] (41) 
the shifted Z-th band of H By 



Aj(e) = (\21-2 - ex, A a -i - ex) (42) 

we denote the shifted l-ih gap. 

As e — > the asymptotics of g(x, E) turns out to be essentially different at 
E G Aj(e) and at E G A/(e). Inside the shifted bands A/(e) 



g(x, E) = (ME)fi(x, E) + Aifi{x, E))[l + 0(e)]. (43) 

Here 

A t (E) = e ir 'Fi_i(E) (44) 

and 

f l (x,E) = (\(k)y 1 / 2 exp 



E+ex 



J (lfc(A) + $ (A)) dX 

^21-2 



<f>(x,\) (45) 



The phase factor rj in (44) can be written down explicitly but we do not need it 
below. The function <p(x, A) in (45) is a periodic part of the Bloch solution ip(x, A) 
normalized by the condition (4). The variables k and A depend on (E+ex) according 
to the relation (40) and $(A) is given by (10). 

Remark. The asymptotics (43) is valid outside the £ 1 / 3 -vicinity of the shifted 
edges A; — ex [14, 16]. In the vicinity of these points one has to use the local 
asymptotics [18]. 



11 



In the shifted gaps Ai(e) outside the e 1 ^ 3 - vicinity of the shifted edges A; — ex the 
solution g has the form 

v, E) = (H^hrix, E) + G^h+ix, E))[l + 0(e)] (46) 



9{x, 



The coefficients Hi(E) and Gi(E) are related with coefficients Ai, Ai as follows 

gMr/^w-a;,)^ + e -in/4 e -i(w-oj t )n ^ 



Hl = ~ 



d = [e^^e^-^Ai + e™ ^ e'^^ A x 

Recall that u and uii are given by (16) and (17) respectively. 
The functions hf(x,E) have the form 



(48) 



hHx,E) = \\(k)\- 1 / 2 exp 



E+ex 



± f (-|9*(A)|+$(A)dA 



[x,X)) (49) 



The function /i + exponentially increases and h exponentially decreases as e — > 
inside the interval Aj(e). 



6 The asymptotics of spectral projections kernels. 

Let us denote by e e (x, y, AE) the kernel of the spectral projection of the operator 
H on the interval AE. Due to (6) it has the form 

1 f 9(x,E)g(y,E) 
e e (x,y,AE) = -J |M(E)|2 (50) 

Theorem 3. Let A' t (e) C A t (e) be a subinterval of the shifted l-th band A t (e) such 
that it does not contain the e 2 ! 3 -vicinity of the edges and A\(e) — > A; as e — > where 
Ai is the l-th band of H . Then at bounded \x\, \y\ 

\ime £ (x,y, A{(e)) = e ; (x,y) 

where ei(x,y) is the kernel of the spectral projection of the operator H on the l-th 
band A; = A/(0). 

Proof. Let us rewrite (45) in the form 

f l (x,E) = (\(k))- 1 ^(x,\y s ^ (51) 

where 

K 

S(x,E,e) = -(E + \ 2l „ 2 )k- J (-A(fc) + dfc. 
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Since k depends on (E + ex) accordingly to (40) then on the basis of (43), (51) one 
has that at E e Ai(e) 

g(x, E)g(y, E) = (T, + T 2 )(l + 0(e)) (52) 
where 

Ti = m^mm lMEW ( 5 3) 

X(k) 

T 2 = 2A" 1 K [^(x, \)^(y, \)e 2iS{x ' E ^A 2 (E)] . (54) 

Consider the contribution of each term in (52) into the kernel e e (x, y, Ai(e)) given 
by (6). Taking into account that \Ai\ = \M\ = \E^\ and using the statement 

of Lemma 3, namely the estimate (38), one obtains 



2,J e> \M(E) ? W ( m 

Make use of (34) and note that A\(e) — > A/ and £ — > A as e — > 0. Then it is 
obvious that the limit of r. h. s. of (55) is exactly the kernel of spectral projection 
ei(x,y) of H on the Z-th band given by (3). Consequently the statement of the 
theorem is valid if 

1 f T 2 (x,E) 

' — >0 as e^O. 



2tt J \M(E)\ 

However this fact is quite obvious due to the rapid oscillations of the phase factor 
exp [iS(x, E, e)] as e — > 0. □ 

Theorem 4. Let A[(e) C A/(e) 6e a subinterval of the shifted l-th gap A;(e) given 
fry (^,2) suc/i i/ia£ i£ does not contain the e 2 ^ -vicinity of the edges and A[(e) — > A; 
as e — > where Ai is the l-th gap of H Q . Then at bounded \x\, \y\ 



hme e (x,y, A{(e)) = 0. 



Proof. At E e Ai(e) the solution g(x,E) is given by (46). For the sake of 
simplicity let us assume that x = y (note that in studying the integrated density of 
states we need namely this case). Then 

g 2 (x,E) = (G 2 (ht) 2 + 2H l G l htK + H 2 (hi) 2 )[l + 0(e)} (56) 

Let us subsequently consider the contribution to e e (x, x, A[) of each term in (56). 
On use of (50) and (56) the contribution of the first term can be written as 

Gi(m<. x ,Ey_ ( _ i{h _- h)+m+m] 



I 



\M(EW 

Me) 
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where 



and 



Here 



and 



h= [ e «("-"«) e 27(*,fl,s) F J-)^ q ( x ,E,e)dE 
J \F N (E)\ 



Me) 



,2j(x,E,e) 



E,e)\Fl-l{Ef 



A 2 (e) 



\F N (E)\< 



E+ex 



-q(x, E, e)dE 



l{x,E,e) = - j \^k{\)\d\ > 

^21-1 



, _. x <p 2 (x,k) 
q{x,E,e) = . exp 
\X(k)\ 

Now make simple estimates of J 2 . Recall that 



E+ex 



2 J $(A)dA 

-^2i-l 



F N {E) 



+ 0{a)). 



^21 



Then notice that at E e Aj(e): 

j(x,E,e)<- J \$tk(\)\d\ = 

end hence 



3 [2 7 (x,£;,e)] < 2tt W 



1 



1 - e- 4CT < 

Inserting (57), (58) into the expression for I 2 one obtains 

HE) 



i 



<3 



tanh a i 



J (1 — e 4<Ti ) sin 2 (c<j — c^) + tanh 2 a"; cos(c<j — uoi) 



dE, 



(57) 



(58) 



where Q > \q(x,E,e)\. Expanding the latter expression in powers of o\ as 07 — > 
one gets 



< - 



4 7 sin 2 (c<j — cj;) 



(l + 0((7,)) = 



14 



eaQ 

47T 



/ A 2 /7T \ / A 2 /-l7T 
Ctg UJl - Ctg | ■ LOl 



ea 



ea 



Thus 



J 2 — > as £ — > 0. 
The integral Ii also vanishes as £ — > 0. Indeed, due the chain of estimates 



IMP 



JV 



= (l + 0(,)) = 



JV 



;i + 0(e))= -5,(^(1 + 0(6)) 



it differs from the integral J 2 by the presence of the rapidly oscillating factor. 
Thus we have shown that the contribution of the first term in (56) into the ker- 
nel e £ (x, x, E) vanishes as e — > 0. 

Consider briefly the contribution of the second term in (56). Notice that 



and hence 



Hid = - (e^-^A? + 



\M(E)\ 2 dE 



-2i(u—u>i 



A«(e) 



because each term in this integral contains the rapidly oscillating factor which causes 
the vanishing of the limits as e — > 0. 

Finally, the contribution of the third term in (56) into spectral projections on 
Aj(e) also tends to 0. It becomes quite clear if one compare this term with the first 
one. Instead of exponentially growing factor e 2l ^ x,E,e ^ the third term exponentially 
decreases as e~ 2 ^ x,E,e \ Since we have proved that contribution of the first term 
vanishes the vanishing of the third one becomes obvious. 



7 Density of States and Integrated Density of States 

Along with the integrated density of states (5) one can consider the density if states 

dN £ {E) 



n e {E) 



dE 



In our case it has the form 



n 



where 



a 

■(E) = - f/3(x,E)dx, 
a J 
o 

1 g 2 (x,E) 



P(x,E) = 



2tt \M(E)\ 2 ' 
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The density of states has no regular limit as e — > 0. Indeed as it has been shown 
in Sec. 4 the spectral density \M(E)\~ 2 is a sort of periodic 5-type sequence and its 
limit is defined only in a weak sense (31). The product g 2 (x,E) due to (52), (54) 
also has no regular limit. 

However the limit of the integrated density of states (5) as e — > is well defined. 
This is the direct consequence of the results of the previous section. Indeed if 
E e A' m+1 (e) (for the notations see theorem 3) the integral (6) at y — x can be 
presented in the form 

e e {x,x,E) = Y, / (3{x,E)dE + ]T / (3{x,E)dE+ 



E 2m 



J P(x,E)dE + J2 J P(x,E)dE. (59) 

E' j=0 Vf 



Here A^(e) and A^(e) are subintervals of shifted bands and gaps described in theo- 
rems 3 and 4 respectively. By we denote the e 2//3 -vicinity of the points Xj — ex 
and E' 2m > \ 2m — ex E V 2 £ m . 

As it follows from theorem 3 at the zero bias limit 

m ,. m 

£ J (3(x,E)dE - 'Eejfax), (60) 
where ej(x,x) is given by (3) and 



jp( X , E)dE (6i) 

E2m A2m 

The consequences of theorem 4 is the limit 

m . 

/ P(x,E)dE^>0. (62) 



A'( £ ) 



One also can prove that 



/ 



(3{x ) E)dE -> 0. (63) 

v? 



To this end one needs a more detailed analysis of the behaviour of (3(x,E) in the 
£ 2 / 3 -vicinity of the shifted edges Xj — ex. This can be done explicitly, however it is 
out the scope of the present paper. 
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Finally inserting the limits (60) (61) (62) (63) into (59) one obtains that at zero 
bias limit 

e £ (x,x,E) - f:ej(x,x) + - f ^ {X \ X ^ X)] dX ^ e(x,x, X). (64) 

\L x ^ 

Here E belongs to the m + 1-th shifted band. 

Analogously one can prove that if E belongs to the shifted m + 1-th gap at the 
zero bias limit 

m 

e £ (x,x,E) -> ^2ej(x,x) (65) 
j'=i 

Since all the limits proved above are valid for arbitrary x, one can integrate them 
w.r.t. x within arbitrary finite interval, say [0, a] 
Thus we came to the following statement. 
Theorem 5 In the zero bias limit e — > 0: E = X and 

N £ (E) -> N(X). 

Here the integrated density of states N £ (E)is given by (5) and N(X) is defined by (2). 
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